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Abstract 

We present a systematic study of twist-3 light-cone distribution amplitudes of K* and 
mesons in QCD. The structure of SU(3)-breaking corrections is studied in detail. Non- 
perturbative input parameters are estimated from QCD sum rules. As a by-product, we 
update the parameters describing the twist-3 distribution amplitudes of the p meson. We 
also review and update predictions for the twist-2 distribution amplitudes of p, K* and (j). 
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1 Introduction 



In recent years SU(3)-symmetry-breaking in processes involving light vector mesons has at- 
tracted increasing interest. For some of these processes, for instance the heavy-meson decays 
B ^ fTf vs. B ^ iir*7, the uncertainty in SU(3) breaking is presently the dominant source of 
theoretical error [HIS]- These and related decays, like B — ^ {p, K*)i~^i~ , are dominated, at 
short distance, by flavour- changing neutral-current transitions which are heavily suppressed 
in the Standard Model (SM) (they occur only at loop-level) and hence are very sensitive to 
potential effects from new physics. All these decays will be studied in detail at the LHC, 
with the aim, in the case of the discovery of new physics in the TeV range, to elucidate its 
flavour structure, or, in the case of continued absence of new particles in direct searches, to 
constrain their possible masses and couplings. In any case, good theoretical control over the 
SM predictions of such decays is vital. The current theoretical approaches to describe them 
all rely, in one way or the other, on their interpretation as hard exclusive reactions, with the 
hard scale set by the heavy quark or meson mass which leads to an expansion in terms of 
the inverse hard scaleQ 

In a perturbative framework, the method of choice for calculating matrix elements of 
B —* light meson transitions is QCD factorisation, which enters QCD sum rules on the light 
cone [4], QCD factorisation for non-leptonic and radiative B decays [S] and perturbative 
QCD factorisation [6]. One important ingredient in these calculations are light-cone hadron 
distribution amplitudes (DAs) which describe the momentum-fraction distribution of partons 
at zero transverse separation in a particular Fock state, with a fixed number of constituents. 
DAs are ordered by increasing twist; the leading-twist-2 meson DA 02; a/, which describes 
the momentum distribution of the valence quarks in the meson M, is related to the meson's 
Bethe-Salpeter wave function ^m.bs by an integral over transverse momenta: 



Here u is the quark momentum fraction, Z2 is the renormalisation factor (in the light-cone 
gauge) for the quark-field operators in the wave function, and fi denotes the renormalisation 
scale. For pseudoscalar mesons, one has one twist-2 DA, whereas for vector mesons there 
are two, and 0|.jvf, one for each independent polarisation state of the vector meson, 
transverse and longitudinal, respectively. In this paper we study the twist-3 distribution 
amplitudes of the vector mesons p, K* and 0, with a particular emphasis on SU(3) (and 
G-parity) breaking effects, and also update earlier results on twist-2 parameters. We do 
not differentiate between p and u mesons as their DAs only differ by the numerical values 
of hadronic parameters which, using the currently available theoretical methods, coincide 
within errors. Our paper is an extension of Ref. |7j to vector-meson DAs and finalises the 
preliminary results for twist-3 parameters quoted earlier in Refs. [U 121 [8] . The results for 
twist-4 DAs will be published elsewhere. 

The study of vector-meson DAs has attracted less attention than that of pseudoscalar 
DAs. The leading twist-2 DAs of p have been investigated in Ref. [9] , correcting a mistake in 
the earlier literature ^L0\. The structure of twist-3 DAs of p, K* and and their relation to 

^This also applies to lattice calculations, see Ref. [Sj. 
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twist-2 DAs have been studied in Ref. |TT], but did not include all SU(3)-breaking effects. In 
this paper, we complete the analysis of Ref. [11] by including all G-parity and SU(3) breaking 
corrections for K* and and also providing numerical values of all hadronic parameters. 

Our paper is organised as follows: in Sec. [2] we introduce notations and review the status 
of the twist-2 parameters. In Sec. [3] we provide parametrisations of all twist-3 DAs to NLO 
in chiral expansion and including all G-parity breaking contributions. In Sec. H] we discuss 
numerical models for these DAs, based on the results from QCD sum rules. We conclude 
and summarise in Sec. [51 The appendices contain a discussion of the renormalisation-scale 
dependence of the twist-3 parameters, which is also affected by SU(3)-breaking corrections, 
and the QCD sum rules used to derive numerical values for all parameters. 



2 General Framework and Twist-2 DAs 

In this section we introduce notations and review the status of twist-2 parameters. 



2.1 Kinematics and Notations 



Light-cone meson DAs are defined in terms of matrix elements of non-local light-ray operators 



extended along a certain light-like direction 
and the meson state. We adopt the generic notation 



0, and sandwiched between the vacuum 



(2.1) 



and 

v[/A^,(«), ... (2.2) 

for two-particle and three-particle DAs, respectively. The superscript A denotes the polar- 
isation of the vector meson: A =|| (_L) for longitudinal (transverse) polarisation. The first 
subscript t = 2, 3, 4 stands for the twist; the second one, M = p, K*, . . ., specifies the meson. 
For definiteness, we will write most expressions for K* mesons, i.e. sq bound states with 
q = u,d. Whenever relevant, we will include quark mass corrections in the form ± m^, 
which allows one to obtain the results for by iris. We do not include the uj, as all 

formulas for DAs coincide with those of the p; the difference is in the numerical values of the 
hadronic parameters which, at least in the framework of QCD sum rules, coincide with those 
for the p except for small differences due to the difference in meson masses. The variable 
u in the definition of two-particle DAs always refers to the momentum fraction carried by 
the quark, u = Ug, whereas u = 1 — u = Uq is the antiquark momentum fraction. The set 
of variables in the three-particle DAs, a = {01,02,0:3} = {os,Og, o^}, corresponds to the 
momentum fractions carried by the quark, antiquark and gluon, respectively. 

To facilitate the light-cone expansion, it is convenient to use light-like vectors and 
instead of the meson's 4-momentum and the coordinate x^: 



I J. Pi J. 



1 



m 



K* 



xP — \ (xP)^ — x'^m 



X, 



2 2 n 
X mj^. 

4:{zpy 



1 

-p,- + Oix' 
2 zp 



" 2 ^ pz 



(2.3) 
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The meson's polarization vector e'-'^-* can be decomposed into projections onto the two light- 
like vectors and the orthogonal plane as follows: 

We also need the projector g^^ onto the directions orthogonal to p and z, 

= 9f^U - — {Pi^Zy + PyZ^) , (2.5) 

pz 

and will often use the notations 

a, = af.z'', bp = h^p^ (2.6) 

for arbitrary Lorentz vectors and 6^. 

The dual gluon field strength tensor is defined as Gfj,y = ^e^ypaG^". Our convention for 
the covariant derivative is = 9^ — igA^. Sometimes, a different convention for the sign of 
g is used in the literature, with = d^ + igA^. The sign of g is relevant for all three-particle 
twist-3 DAs, see Tab. [H 



2.2 Conformal Expansion and the Structure of SU(3)-Breaking 
Corrections 

A convenient tool to study DAs is provided by conformal expansion, see Ref. [I2] for a 
reviewH The underlying idea is similar to partial-wave decomposition in quantum mechanics 
and allows one to separate transverse and longitudinal variables in the Bethe-Salpeter wave- 
function. The dependence on transverse coordinates is formulated as scale dependence of 
the relevant operators and is governed by renormalisation-group equations, the dependence 
on the longitudinal momentum fractions is described in terms of irreducible representations 
of the corresponding symmetry group, the coUinear conformal group SL(2,M). 

To construct the conformal expansion for an arbitrary multi-particle distribution, one first 
has to decompose each constituent field into components with fixed Lorentz-spin projection 
onto the light-cone. Each such component has conformal spin 

J = lil + s), 

where I is the canonical dimension and s the (Lorentz-) spin projection. In particular, 
1 = 3/2 for quarks and / = 2 for gluons. A quark field is decomposed as ip^ = A+ip 
and ip- = A_ip with spin projection operators A+ = '-^p-^z/ i^pz) and A_ = -^zlpl i^pz), 
corresponding to s = +1/2 and s = —1/2, respectively. For the gluon field strength there 
are three possibilities: Gz± corresponds to s = +1, Gp± to s = —1, and both G±± and G^p 
correspond to s = 0. Multi-particle states built of fields with definite Lorentz-spin projection 
can be expanded in irreducible representations of SL(2,M) with increasing conformal spin. 

^See Ref. |13) for an alternative approach not based on conformal expansion. 
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The explicit expression for the DA of an m-particle state with the lowest possible conformal 
spin j = ji + . . . + jm, the so-called asymptotic DA, is given by [1] 



i / \ r(2ji + ■ ■ ■ + 2jm) 2n-l 2i2-l 27„,-l /o ^7^ 

r(2ji) ■ ■■T{2jm) 

Multi-particle irreducible representations with higher spin j + n,n = 1,2, .. ., are given by 
polynomials of m variables (with the constraint YlT=i'^k = 1 ); which are orthogonal over 
the weight function (12 .Tp . For the twist-2 and 3 two-particle DAs these are Gegenbauer 
polynomials, whereas the twist-3 three-particle DAs get expanded in Appell polynomials. 

In this paper we are particularly interested in SU(3)-breaking corrections to DAs. These 
corrections come from different sources: 

• SU(3) breaking of hadronic parameters: these effects are partially known for twist-2 
parameters, see Refs. [IDl [HI [15], but have not been studied for twist-3 parameters 
before; 

• G-parity breaking parameters: these are of parametric order — mg and vanish in 
the limit of equal quark mass, i.e. for p and 0. For twist-2 DAs, these have been 
calculated, to lowest order in the conformal expansion in Refs. [IHl [Til [151 [161 [17] ; they 
are unknown for twist-3 DAs@ 

• explicit quark mass corrections in m^img to DAs and evolution equations: these affect 
only higher-twist DAs and are induced by the QCD equations of motion (EOM) which 
relate twist-3 DAs to each other and to twist-2 DAs, see Sec. [3] The mass corrections 
to vector meson DAs have been calculated to twist-3 accuracy in Ref. [H]; the effect 
on the evolution of DAs under a change of the renormalisation scale so far has only 
been investigated for pseudoscalar DAs [7]. 

We shall study all these effects in this paper. 

Let us now see how these corrections affect twist-2 DAs. 

2.3 Twist-2 Distributions 

The twist-2 DAs of K* mesons are defined in terms of the following matrix elements 

of non-local operators = 2u — 1) [llj: 



{0\qixM-x)\K*{P,X)) = flniK* \-^P,j^ 



1 e^^^x 



^The results given in Refs. [H [2l [8] were preliminary versions of those obtained in this paper. 
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{0\q{x)a^,s{-x)\K*{P,X)) = 
1 



1 
4 



e^-^^x 2 



' II 1 1 



Px 



Jo 



[2.9) 



All other DAs in the above relations are of twist 3 or 4. We have neglected all terms in the 
light-cone expansion which are of twist 5 or higher. The normalisation of all DAs is given 
by 

(2.10) 



du (j){u) = 1 . 

The above DAs are related to those defined in Refs. [HI [18] by 



h-K' — y± 



(2.11) 



The conformal expansion of 02 reads 

f oo 



(2.12) 



n=l 



in terms of the (non-perturbative) Gegenbauer moments ai' and the Gegenbauer polyno- 
mials Cn ■ To leading-logarithmic accuracy, the renormalise multiplicatively as 

af(/.^) = L-^"/(^^o)a„(/.g), (2.13) 

where L = as(/i^)/as(/io), = (33 — 2Nf)/3, and the anomalous dimensions 'jn^ are given 

by m 



7p) = 8CFi^ijin + 2)+^E-- 
7„^(°) = 8Cf (^|J{n + 2) + 7i. - 



4 2(n + l)(n + 2) 
3 



To next-to-leading order accuracy, the scale dependence of the Gegenbauer moments is more 
complicated and reads ^20j 



n-2 



(2.14) 



fc=0 



where 



1 Svr/?^ 



with /5i = 102 — (38/3) A^/; 7!^'' are the diagonal two- loop anomalous dimensions, which 
have been calculated, for the vector current, in Ref. [21], and, for the tensor current, in 
Ref. [22]. The mixing coefficients d^^l, k < n — 2, are given, in closed form in Ref. [20], for 
the axial vector current; the formulas are valid for arbitrary currents upon substitution of 
the corresponding one-loop anomalous dimension. 

For the lowest moments n = 0, 1, 2 one has, explicitly: 

lim iim 23110 512 ,^ nm 34072 830 

0^'"^124-8A^,. ,i'^'.'^-^N„ (2.15) 

and 

11(1) _ 35 20 - 3/9o . ;-50/{9/3o)-n 

j.(i) ^ 28 16-3^ . _ ^4o/(9/3o)-n , (2.I6) 
^° 9 40 - 9/3o ^ ^ 

Let us now review the numerical values of the twist-2 parameters, to NLO in conformal 
spin. The longitudinal decay constants of the charged mesons p^, K*^ can be extracted 
from the branching ratios of V~Ut-, whereas fp^^ follow from e+e^ — > V^. A critical 

discussion of the results, including effects of p-u and uj-(f) mixing, was given in Ref. [2] from 
which we quote the following results: 

41 = (216 ± 3) MeV, £ = (187 ± 5) MeV, /i = (220 ± 5) MeV, ^ = (215 ± 5) MeV . 

(2.17) 

The transverse decay constants /-*", on the other hand, cannot be determined from ex- 
periment, but have to be calculated using non-perturbative methods. Currently available 
results include QCD sum rule determinations [9], HH |15l [17] and lattice calculations [23] . 
The corresponding results have been critically reviewed and averaged in Ref. [2], with the 
following results: 

/^^ = (165 ±9) MeV, = (151 ±9) MeV, = (185 ±10) MeV, = (186±9) MeV . 

(2.18) 

Let us now turn to the Gegenbauer moments a|'^. At present, there are no lattice 
determinations for any of those, so all available determinations come from QCD sum rules 
[m HSl [m ITTI [2I] or quark models [25]. For mesons with definite G parity (equal mass 
quarks), i.e. p and (j) in our case, af'"*" = 0. For al''^{K*), the results from QCD sum rule 
calculations converge to [17] 

a|(if*)'^=^^°^ = 0.03 ± 0.02 , al{K*Y=^^''^ = 0.02 ± 0.02, 
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As for 02, a2 (p) have been determined in Ref. [9J and reinvestigated,recently, in Ref. [Tj, 
using the updated hadronic input collected in Tab. |Al The resulting values 

^|^^^^=iGcV = 0.15 ±0.07, a|(p)''=2^<=^ = 0.10 ±0.05, 

4(P)''^'^°^ = 0.14 ± 0.06 , a^(p)^=2G<=v = 0.11 ± 0.05 , (2.20) 

are slightly smaller than those quoted in Ref. [5]. The value of a\''^{K*) has been determined 
in Ref. pTl [T5] and reinvestigated in Ref. [Ij. The result is 

aj(fs:*)^=^^°^ = 0.11 ±0.09, a|(ir*)^=2^'^ = 0.08 ±0.06, 

a^(ir*)^=^^^^ = 0.10 ±0.08, a^(ir*)^=2^'=^ = 0.08 ±0.06. (2.21) 

The corresponding parameters of the have received far less attention: Ref. [TT] quotes 
a|'^(0)^=^^^^ = 0±0.1 and Ref. [2] 0^(0)^=^^*^^ = 0.2 ±0.2. In this paper, we evaluate the 
sum rules collected in App. [Aj which include all relevant corrections in and in particular 
the radiative corrections to the quark condensate term in ms{ss), and find 

4(0)'^=^^"^ = 0.18 ± 0.08 , a|(0)''=2^"^ = 0.13 ± 0.06 , 

0^(0)'^=^^'^^ = 0.14 ±0.07, a^(0)^=2G"V = 0.11 ±0.05. (2.22) 

In summary, it is probably fair to say that all known determinations of 0^2" point at 

fairly small values at 1 GeV and that within the present accuracy a^^) = '^1(2)- 
Let us now turn to twist-3 DAs. 



3 Twist-3 Distributions 

To twist-3 accuracy, there is a total of four two-particle DAs and three three-particle DAs 
whose mutual interrelations have been unravelled in Ref. [IT], including quark mass correc- 
tions. The crucial point in constructing higher-twist DAs is the necessity to satisfy the QCD 
EOM which yield relations between physical effects of different origin: for example, using 
EOM, the contributions of orbital angular momentum in the valence component of the wave 
function can be expressed in terms of contributions of higher Fock states. An appropriate 
framework for implementing these constraints was developed in Ref. [14] : it is based on the 
derivation of EOM relations for non-local light-ray operators [26j, which are solved order by 
order in the conformal expansion; see Ref. [T2| for a review and further references. In this 
way one can construct self-consistent approximations for the DAs, which involve a minimum 
number of hadronic parameters. The EOM relations relating twist-2 and -3 DAs of vector 
mesons were derived in Ref. [TT], including all quark-mass corrections. What is new in the 
present paper is the inclusion of G-parity breaking corrections to three-particle DAs, which, 
via the EOM relations, also impact on the two-particle DAs. Based on these relations, we 
derive, in this section, complete formulas for all twist-3 DAs to NLO in the conformal expan- 
sion, including all G-parity breaking effects. A non-zero quark mass also induces a mixing 



7 



of twist-2 parameters into those of twist-3 under a change of the renormahsation scale. We 
also derive the corresponding scaling relations. 

Let us start by defining the relevant DAs. The two-particle twist-3 DAs (p^.'x* have 

already been defined in Eqs. (12. 8p and (12. 9p . There are two more two-particle DAs, ip^.x*, 
defined a^ 

{0\q{zh,^,s{-z)\K*{P,X)) = i/im^.e/-M'W/5^'^«e^^^^V^3V*(«), (3.1) 

{0\q{z)si-z)\K*iP,\)) = -z/i.(eWz)m^. [' due'^^^4.^^,{u) . (3.2) 

Jo 

The normalisation is given by 



du ^J] y, {u) = l- -j^ — , (3.3) 

fK* 



which differs from Ref. [H], where all DAs were normalised to 1. The reason is that in 
[TT] we implicitly expanded the normalisation factor 1/{1 — (/J[-l"'"V/x*"^)("^s + ^ij)/"^x*} in 
powers of + m^, whereas in this paper we keep the full dependence on the quark masses. 
There are also three three-particle DAs of twist 3: 

(0|g(%G^,(t;z)7.755(-^)|/r(P,A)) = fl.mK^{pzfe%^.^J,{v,pz) + ..., 
{mz)9Gp.{vz)t^A-^)\K*{P. A)) = fl,mK^{pzfef^^lj,{v,pz) + . . . , 

{mz)gG-Avz)a-^ps{~z)\K*{PA)) = fk^ml,{e^'^ z){pz)^i^K.{v,pz) , (3.4) 
where the dots denote terms of higher twist and we use the short-hand notation 

^(„.p.)^/l>^e->~)^fe) (3.5) 



with J^{a) being a three-particle DA. a is the set of parton momentum fractions a 
{«!, a2, as} and the integration measure Va is defined as 



J'Da = J daida2daz 5{1 — aj) . 



(3.6) 



As discussed in Ref. [H], all these DAs are interconnected by the QCD EOM. The analysis 
of these EOM and the resulting relations including quark mass corrections in iris ± f^iq is the 
subject of Ref. [11], so we do not repeat it here, but just quote the results: 



1,„. . /-^ 1 



ipl-K^i^) = u I dv-T{v)+u / dv-T{v), 



^In the notations of Ref. ^i;K' ^ i^-ifk- / fl')ims+mg)/mK'}g'f \ 4'1,k' ^ {^-(fl-/ fk')ims + 



dv - T(v) 



[ dv-T(v) 

Ju V 

/ dai I da2— ^^.k'{ql) 
du Jo Jo as 



f^. rus + mq II 



(3.7) 



with 



fl, TTLs + m„ 



2^du 



1/1. rris - rriq d j 

[u]--— -r n:K*W 



2 /k* riiK* du 



"3 



dai 



da-2 



(3.8) 



and 



1 /"^ 1 

'ijj^^j^*{u) = u I dv-Q,{v)+u / dv-Q.{v) 

'o ^ Ju 



dv-Vtiv) + / dv-Vtiv) 

^ in 



+ 



K* 



niK* 



+ 



du 



/ dai I da2 — ^1-k*{(^) 

Jo Jo Oi3 



+ I dai J da2 — ( -r— + -7— 1 ^l-h-A^) 



03 \dai ' da2j 



(3.9) 



with 



Vt{u) = 20^^. M + ^ ^ — 



fl, du 



d r f"^ 1 f d d 

+ 2— / dai I da2 — "i ^ h 0^2 ^ — \^\ k*{^ 

du Jq Jq as \ dai da2 



d r r I f d d 

+ 2 — / dai / da2 — ( ai 0:2 - — ) $3.7^. (a) 



du 



"3 



dai 



da-) 



(3.10) 



The twist-3 three-particle DAs correspond to the hght-cone projection 'jzGz± and <7j_zGj_z, 
respectively, which picks up the s = | component of the quark fields and the s = 1 component 
of the gluonic field strength tensor. According to (12.71) . the (normalised) asymptotic DA is 
then given by 3Q0aia2a^. To NLO in the conformal expansion, each three-particle twist-3 
DA involves three hadronic parameters, which we label in the following way: (, k are LO and 
uj, A NLO parameters. ( and u are G-parity conserving, whereas k and A violate G-parity 
and hence vanish for mesons with quarks of equal mass, i.e. p and (p. We then have 

^i,K*i^) = 360aia2a3 |i* + " "2) + ^Ik-^ (7tt3 - 3)| , 



"^ix^ioi) = 360aia2a3 jcL* + A!l^*(ai - 02) + c^L*^ (''"s - 3)| , 

H;K*i^) = SeOaiasa^ j^J-^^ ^ - ^2) + A^^,^ (Tag - 3)| . (3.11) 

The relation to the parameters used in Ref. \TT\ is = C|, C3 = ^^3/14, (J = 00^ /lA, 
Clojf'Q = ujI', G-parity breaking terms were not considered in Ref. |Tr|. For equal mass 

quarks, ^^.'J^* are antisymmetric under ai ^ 02, whereas $3.^, is symmetric. 

All these parameters can be defined in terms of matrix elements of local twist-3 operators. 
For chiral-odd operators, for instance, one has 

{0\qa,^gG,^s\K''{P,X)) = f^.ml,{e^^h){pz)KiK. , 

{0\qtD.a,i^gG,^s - qa,i:gG,^iD.s\K*{P,X)) = f^.ml,{e^^h){pzf ^iuij,. ; (3.12) 

the formulas for chiral-even operators are analogous. Numerical values for these parameters 
can be obtained from QCD sum rules and will be discussed in Section HI 

Of the parameters in (13.111) . CJk*; '^s/j'*; ^sk* ^^nd A^^. renormalise multiplicatively in 
the chiral limit, the others mix with each other. For non-zero strange quark mass, there 
is additional mixing with twist-2 parameters. Here, we write down explicitly only the RG- 
improved relations for the above 5 parameters; a full discussion, including also A|, c^I, is 
given in App. |Al The relations can be written in compact form as 

3 

i=i 

with the LO scaling factor L = as{fJ''^)/c(s{f^o)- If there is a fiavour threshold /ith between /iq 
and /i changing the number of active fiavours from Uf to Uf + 1, then one has to replace 

lV/^o ^ («(/i^)/«(/it\))^/^"^"^+^H«(/x,\)/a(/.^))^/^°("^-). 
The parameters in (I3.13P are given by: 

P = {/I-'CIe-*) /l'*'^3i<:*' fK'^^SK'y fK*^3K*J /if'Ag;^,}, 

Qi(2) = — ^ {rris ± TUg, {rus =F mg)a^, {rUs ± mq)a^} , 
rriK* 

/" - II II 

Q3,5 = — ^ {ms - rUg, {nis + mg)a\, {nis - mg)al} , 

/" • II II 

Qi = — ^ {rris + rrig, {rris - mg)al, {rUs + mg)al} , 
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77 77 55 73 104 



7p = i TT' TT' TT' TT' 



9 9 9 9 9 



(7q) 



1,2 



16 88 



(7(9)3,4,5 



68 86 



2 

29 
2 

29 
4 

r9 

14 

37 
1 

V~85 



6 

25 
6 

'25 
12 

65 
42 

'25 
1 

~5 






^ 

12 

13 
4 

15/ 



(3.14) 



Implicit formulas for the remaining 4 parameters in (13. lip can be found in App. \M Numerical 
values will be given in the next section. 

Using (13. lip , and the corresponding relations for twist-2 DAs, one obtains expressions 
for the twist-3 two-particle DAs, which are valid to NLO in the conformal expansion. As 
discussed in Ref. [llj, the structure of this expansion is complicated by the fact that these 
DAs do not correspond to a fixed projection of the quark fields' Lorentz-spin s. The resulting 
expansion is in Cn^'^iO "^t-K* ^'^d Cn^'^iO for (f)'^.'^*- 

+ {j - I - 3) + ^ u^ij,, (3 - 30e + 35e^) 

+ - ^'^^1 ^ |l + 8^al + 3(7 - 30uu)al + ^lnM(l + 3a} + 6a|) 
2 Ik* ^ 



— ^lnM(l — 3af + 6a|)| 



3 TUs - TUg fj^, 
2 rriK- fk* 



^{2 + 9^a! + 2(11 - 30nn)a5 + lnM(l + 3a! + 



+ In u{l — 3af + 6a|) | 



(3.15) 
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+ 3 ^'^^1 Ij<1 Juu(l + 2{a| + 3(7 - 5uu)al) +u\nu(l + 3a} + 64) 
TTiK- fx* y 

+u In «(! — 3af + 6a|) | 
-3 ^ {uu{'^a\ + 10^a|) + u \nu{l + 3a} + 6a|) 

— M lnM(l — 3a} + 6a|)| , (3.16) 
i'i;K*{u) = 6uu |l + Qa} + yi*) Cl^\0 

+ 6 ' " ^ {uu{2 + S^ai- + 2(11 - 10MM)a^) +ulnu{l + 3a^ + 6a^) 

+u \nu{l — Sa^ + 6a^)} 
-6 ij<l [uu{9a^ + lO^a^) + u \nu{l + 3a^ + 6a^) 



-ulnM(l -3a^ + 6a^)} , (3.17) 



+ + l + {^-al + 5CL.| (3^^ - 1) + js^. - 1^ 



+ - !!^i±I^ 4^ {2 + g^a^ + 2(11 - 30nn)a^ 

+ (1 - 3a^ + 6a^) \nu + (1 + 3aj- + Qa^) \nu} 

2 m^. /II , 

+ (1 + 3a^ + 6a^) In u - (1 - 3a^ + 6a^) Inu] . (3.18) 

These expressions are our final results for the two-particle twist-3 DAs and supersede those 
given in Ref. [TI] where G-parity violating terms in and A3 were not included. The 
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P 


K* 







fi = 1 GeV 


H = 2 GeV 


fi = 1 GeV 


/i = 2 GeV 


/i = 1 GeV 


/i = 2 GeV 


a\ 








0.03(2) 


0.02(2) 
















0.04(3) 


0.03(3) 








al 


0.15(7) 


0.10(5) 


0.11(9) 


0.08(6) 


0.18(8) 


0.13(6) 




0.14(6) 


0.11(5) 


0.10(8) 


0.08(6) 


0.14(7) 


0.11(5) 




0.030(10) 


0.020(9) 


0.023(8) 


0.015(6) 


0.024(8) 


0.017(6) 










0.035(15) 


0.017(8) 








nil 


-0.09(3) 


-0.04(2) 


-0.07(3) 


-0.03(2) 


-0.045(15) 


-0.022(8) 










0.000(1) 


-0.001(2) 










0.15(5) 


0.09(3) 


0.10(4) 


0.06(3) 


0.09(3) 


0.06(2) 


a" 








-0.008(4) 


-0.004(2) 
















0.003(3) 


-0.001(2) 










0.55(25) 


0.37(19) 


0.3(1) 


0.2(1) 


0.20(8) 


0.15(7) 










-0.025(20) 


-0.015(10) 









Table 1: Twist-2 and -3 hadronic parameters at the scale = 1 GeV and scaled to up ^ = 2 GeV, 
using the evolution equations ()3.13p . The sign of the twist-3 parameters corresponds to the sign 
convention for the strong coupling defined by the covariant derivative = — igA'^t'^; they 
change sign if g is fixed by = 5^ + igA'^t'^. 



terms multiplying \nu and \nu are the first three terms in the conformal expansion of 
d/{du)(j)2.j^,{0) and d/{du)(j)l'.'l^,{l), respectively. Numerical values for the hadronic param- 
eters are given in Tab. [TJ 

4 Models for Distribution Amplitudes 

In this section we compile the numerical estimates of all necessary parameters and present 
explicit models of the twist-3 two-particle distribution amplitudes introduced in the last sec- 
tion. The important point is that these DAs are related to three-particle ones by exact QCD 
equations of motion and have to be used together; this guarantees the consistency of the 
approximation. Our approximation thus introduces a minimum number of non-perturbative 
parameters, which are defined as matrix elements of certain local operators between the vac- 
uum and the meson state, and which we estimate using QCD sum rules. More sophisticated 
models can be constructed in a systematic way by adding contributions of higher conformal 
partial waves when estimates of the relevant non-perturbative matrix elements will become 
available. 

Our approach involves the implicit assumption that the conformal partial wave expansion 
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is well convergent. This can be justified rigorously at large scales, since the anomalous 
dimensions of all involved operators increase logarithmically with the conformal spin J, but 
is non-trivial at relatively low scales of order fi ~ (1^2) GeV which we choose as reference 
scale. 

Since orthogonal polynomials of high orders are rapidly oscillating functions, a trun- 
cated expansion in conformal partial waves is, almost necessarily, oscillatory as well. Such 
a behaviour is clearly unphysical, but this does not constitute a real problem since physical 
observables are given by convolution integrals of distribution amplitudes with smooth coef- 
ficient functions. A classical example for this feature is the 77*-meson form factor, which is 
governed by the quantity 



where the coefficients a, are exactly the "reduced matrix elements" in the conformal expan- 
sion. The oscillating terms are averaged over and strongly suppressed. Stated otherwise: 
models of distribution amplitudes should generally be understood as distributions (in the 
mathematical sense). 

We give all relevant numerical input parameters for our model DAs in Table [H at the 
scale /i = 1 GeV, which is appropriate for QCD sum-rule results, and, using the LO and 
NLO scaling relations given in Sees. 12.31 and [3l at the scale /i = 2 GeV, in order to facilitate 
the comparison with future lattice determinations of these quantities. The mixing of K* 
and parameters with operators of lower twist depending on is numerically small. In 
evaluating the sum rules, we have chosen the values of the continuum threshold sq as given 
in App. O The sum rules are actually rather insensitive to that parameter, due to the 
smallness of the perturbative contribution, but are not very stable in the Borel parameter 
M^, which is the reason for the large uncertainties in Tab. [H The biggest contribution to 
G-conserving parameters comes from the gluon condensate. For G-breaking parameters, on 
the other hand, this contribution is suppressed by a factor and, as a consequence, all 
G-breaking parameters are considerably smaller than the G-conserving ones. The table also 
shows that SU(3) breaking is relevant for all parameters. 

At this point we would like to compare our results with those available in the literature. 
We have discussed the twist-2 parameters already in Sec. 2.3. As for twist-3, only results for 
p are available. A long time ago, the chiral-even parameters were determined, in Ref. \27\ . 
as 



A comparison with Tab. [T] shows that these values agree quite well with ours, although we 
think that the uncertainty of (Ip was underestimated in |27j. A value for was obtained 
in Ref. [H]: a;^(l GeV) = 0.3 ± 0.3, which is a bit smaller than our result. 

In Fig. [T] we plot the longitudinal twist-3 two-particle DAs 0| and ipl for the p meson, 
assuming massless quarks, and for the K* and (j) mesons, together with the corresponding 
asymptotic DAs. Figure [2] shows the transverse DAs (p^ and ip^. The figures show that 
quark-mass corrections significantly modify the end-point behaviour of ^f'"*", where they 
induce a logarithmic end-point divergency, even if the contributions of gluonic operators are 
neglected. This is not a problem because, as mentioned above, the DAs themselves need 
not be finite, it is only their convolution with perturbative scattering amplitudes that is 




Ci(l GeV) = 0.033 ±0.003, cj" (1 GeV) = 0.2 , cD" (1 GeV) = -0.1 . 



(4.1) 
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0.2 0.4 0.6 0.8 1. 0.2 0.4 0.6 0.8 1. 

U U 

Figure 2: Left panel: cf)^ as a function of u for the central values of hadronic parameters, for 
/i = 1 GeV. Red line: 4>3;p^ green: (I)3-k*i blue: 0^^- Right panel: same for ip^. 



meaningful. The figures also show that the effect of SU(3) breaking (the difference between 
the red and the other curves) is quite pronounced for all DAs, whereas the G-parity breaking 
terms (the asymmetry of the green curve) have only minor impact, which is due to the 
numerical smallness of the corresponding hadronic parameters. 

5 Summary and Conclusions 

In this paper we have studied the twist-3 two- and three-particle distribution amplitudes 
of K* and (p mesons in QCD and expressed them in a model-independent way by a min- 
imal number of non-perturbative parameters. The work presented here is an extension of 
Refs. [TJ [11] and completes the analysis of SU(3)-breaking in vector meson distribution am- 
plitudes to twist-3 accuracy. Our approach consists of two components. One is the use of the 
QCD equations of motion, which allow dynamically dependent DAs to be expressed in terms 
of independent ones. The other ingredient is conformal expansion, which makes it possible to 
separate transverse and longitudinal variables in the wave functions, the former ones being 
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governed by renormalisation-group equations, the latter ones being described in terms of irre- 
ducible representations of the corresponding symmetry group. We have derived expressions 
for all twist-3 two- and three-particle distribution amplitudes to next-to-leading order in the 
conformal expansion, including both chiral corrections 0{ms + mq) and G-parity-breaking 
corrections 0{ms — rriq)] the corresponding formulas are given in Sec. [3l 

We have also done a complete reanalysis of the numerical values of the relevant twist-3 
hadronic parameters from QCD sum rules. Our sum rules can be compared, in the chiral 
limit, with existing calculations for the p |Tll[27]. We have also studied the scale-dependence 
of all parameters to leading- logarithmic, or, if possible, next-to-leading-logarithmic accuracy, 
taking into account the mixing with operators depending on the strange-quark mass m^. Our 
final numerical results, at the scales 1 and 2 GeV, are collected in Tab{T] 

Preliminary versions of our results have already been applied in studies of B ^ (p, K*)'y 
decays [Il[2l[8]. These processes are also sensitive to twist-4 distribution amplitudes which 
we will study in a separate publication. While the parametrisations given in Sec. [3] are 
general, the actual numerical values for hadronic parameters given in Sec. 12.31 and Tab. [T] are 
obtained within the framework of QCD sum rules. We are looking forward to a confirmation 
of these values from independent non-perturbative methods, for instance lattice QCD. 
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Appendices 

A Scale-Dependence of Twist-3 Parameters 

In this appendix we derive the renormalisation-group improved expression fl3.13p for twist-3 
parameters and also give an implicit relation for the scaling of the chiral-even parameters 
^^3, '^s, As and A3. 

Let us introduce the following notation for the relevant quark-quark-gluon operators, see 
Eq. (EUD: 

0^{z,vz,-z) = q{z)a^ugGzu{vz)s{-z) , 

Os{z, vz, -z) = q{z)gG±^{vz)i-f;,s{-z), O^iz, vz, -z) = qiz)gG±^{vz)'j^'y5s{-z) . (A.l) 

In principle it is possible to establish evolution equations for these non-local operators using 
the light-ray operator technique developed in Ref . [25] . The scale-dependence of the param- 
eters in (13. lip then follows from a projection of the evolution equation on the corresponding 
conformal wave. Another approach is to make use of the results derived in the literature 
for the anomalous dimensions of moments of the corresponding nucleon structure functions. 
The explicit relations between these anomalous dimensions, given in Ref. [28], and those 
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needed for our twist-3 parameters have been established in Ref . [TT] . Neglecting quark mass 
corrections, one has, in the notations of Ref. fill. 



o) ; 



II 

2 II 



3X* 



3K* / 



/2 II 



a;. 



Mo 



- ^3^* + ^3^, j 



fx* 



3X» 



+ A 



3K* 



/Has 



+ A 



Mo 



3 



3/<* / 



(A.2) 



where L is the leading- log scaling factor L = as(/i^)/as(/ig). The factor 2/3 comes from the 
relative factors in (13. 121) . The anomalous dimensions are given by 



J- 2 



7 

Ca + 3 Cp , 



Ca + — Cp , 





— C/i -l- — Cp . 
3 6 



3C/ 



2 5 



— C A H — C'f 
3 3 



2 2 



ACa + Icp 
o 

2 2 

— C/i H — Cp 

3 3 



(A.3) 



For massive quarks, the scaling relations receive corrections in ± m^, depending on 
the G-parity of the parameter. These corrections are induced by mixing of the operators in 
(lA.ip with twist-2 operators and can be calculated using the light-ray-operator technique of 
Ref. [26j resulting in the following compact expressions: 



Ol{z, vz, QY" = Ol{z, vz, 0)'^o + , ^ In ^ 



dt 



02(z, vz) — 2t02{z, tvz) 



m„ 



02ivz, 0) - 2t02((l - (1 - v)t)z, 0) 



+ . . . 
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In 4 t dt\'^ 

V 



1-v 



47r Jo 
OUvz, 0)-2tOU{l-{l-v)t)z,0) 



Ol{z,vz) - 2tOl{z,tvz) 



2 /•! 



dt 



V 



O2 (-2, VZ) -2tOl{z,tvz) 



rrir. 



1-v 



Olivz, 0)-2tO^{{l-{l-v)t)z,0) 



+ ... 



(A.4) 



The twist-2 operators in the above relations are given by 

02{az,bz) = q{az)'^ zs{hz) , Oj(a2;, 62;) = q{az)a ^_zs{hz) (A. 5) 

and the dots stand for 0{as) contributions from the twist-3 operators. 

Taking (1A.2I1 and (1A.4P together, one finds that the G-even parameters mix with 

whereas the G-odd ones mix with 

/^l"^{(ms - ruq), {nis + mg)af'^^^\ {nis - mg)a^^"^}. 

For the _L parameters, and (I and k|, the combination of ( 1A.2I) and the projection of ( 1A.4P 
onto the corresponding partial waves results in the scaling relations given in fl3.13p . For the 
remaining parameters, one finds the following relations: 



Mo 



A ~\\ 

rris + niq . 

Ik* 

rriK' 
iTLs — rriq j_ I 

rriK' 

\ — Jk*(^2 

fx* ^Ik* 



J 



fK*^3K* 



rrts + mq j_ 



rriK* 
nis - ruq j_ 



rriK* 



Jk*0'2 ) 



J 



nis - ni„ 



niK' 



^^^^^ fk^aUK*) 
\ rriK' 



\ niK* 

( f" a" X'^' 

fx* ^3K* 



ms - niq 
niK' 



K* 



niK' 
\ niK' 



(A.6) 
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The anomalous dimension matrices are given by 



/ 179 


7 


7 


7 


_ ' 




/ 77 


1 


2 


2 


4 \ 


18 


4 


~9 


15 


3 




6 


3 


45 


5 


15 


1 


77 


2 


2 


4 




7 


179 


7 


7 


2 


3 


y 


~45 


~5 


^15 




4 


Is" 


9 


^15 


~3 








16 

y 
















16 

y 

















8 
















8 





V 











88 















88 

y/ 



The resulting explicit expressions for the twist-3 parameters are rather bulky, so we do not 
give them explicitly. 



B Sum Rules for Twist-2 Matrix Elements 



In this appendix we list and evaluate the QCD sum rules for twist-2 matrix elements of 
the K*. The sum rules for /|^^, including SU(3)-breaking corrections, were calculated in 
Refs. dZlEn], those for a\'^{K*) in Ref. [HI, and those for al'^{K*) in Ref. [15], apart from 
the perturbative terms in and the radiative corrections to the quark condensate, which 
are new. 

For the longitudinal parameters, the sum rules read: 



so 



47r2 



2 (s-m2)2(s + 2m2) a,, 



+ 



IT Att"^ 



1 - e-^°/*^' 



+ 



M2 
167ra, 



1 + _ 1^ _^ 1 ^ ms{qq) 



{qq){ss) 



3M2 9 TT 

167rQ; 



9M4 



81M4 
5 



3 7r M2 
{{qq)' + {-ssf) , 



so 



47^2 



2 [s — m 



2\2 



+ 



5mz 



1 + 



/I ^. / So \ , m2 M2 /M2 \ 

l-7. + ln^ + ^e-/- + E.(- 



18M4 \ TT 

5 ms{ss) 
3 M2 

5 mg(ss) 20 ms{qq) 5 ms{sagGs) 80na 



124 8 
yr ^ 9 



3 



+ 



27 TT M2 



+ 



9 



+ 



81M4 



(B.i; 



(B.2) 
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so 



36M2 \ TT 



+ - 



7 ms(ss) f ^ , a 
3 M2 



1 + — 

TT 



27 



184 25 / , M2 M2 ^. / So \\n 



49 tts ^^(gg) 35 ms{sagGs) 2247ras /,_ >2 , /- \2\ 1127rQ;3 „s 
+ — — — — + [{qq) + {ss) ) - {qq){ss). (B.3) 



27 TT M2 18 M4 81M4 
For the transverse parameters, one has: 



so 



s — miy{s + 2m 



47r2 J TT \9 3 iJ?J 



ms{ss) 



1 ,^^2. 



12M2 ' TT 



M2 
1 

3M^ 



as f 22 2 



ms{sagGs) 



3M2 TT V 9 3 

327rQ; 



81M4 



(B.4) 



5m2 ,Q; 



47r2 ^ 



s — mIY 10 ms{sagGs) 



+ — 



9 



9M4 



_ 5 m^(gg) Lj^^ 



3 M2 



TT 



49 4 / , ^,j^2 ^. / So \ 



M2 M2 

— + — 
fJ' So 



, (B.5) 



so 

47r2 



V S5 QOtT^ TT ^ 



7 /a 



54M2 \ TT 



+ 



7 msjss) , ^ 



3 M2 



TT 



206 16 



27 



20 



m) 

{qagGq) 


= (—0.24 ± U.Ulj GeV 

= "^0 km) 


{ss) = 
(sagGs) = 


(1 - O3) {qq) 
(1-6,) {qagGq) 




= (0.012 ± 0.003) GeV^ 






ml 


= (0.8±0.1)GeV^ ^3 


= 0.2 ±0.2, 


65 = 0.2 ±0.2 



m,(2 GeV) = (100 ± 20) MeV < — > m,(l GeV) = (133 ± 27) MeV 
mgifi) =ms{fi)/R, R = 24.6 ±1.2 



a,{mz) = 0.1176 ±0.002 



a,(lGeV) = 0.497 ±0.005 



Table A: Input parameters for sum rules at the renormalisation scale = 1 GeV. The value of 
is obtained from unquenched lattice calculations with nj = 2 flavours as summarised in [30] , which 
agrees with the results from QCD sum rule calculations [31]. ffiq is taken from chiral perturbation 
theory ^2]- ots{mz) is the PDG average 



49 ms{sagGs) lUnas ,2 , /- \2\ 



(B.6) 



For p and 0, one has af'"*" = 0. To obtain the sum rules for /jl'"*" and a|"'"(0), one has to 
replace the perturbative contributions to the above sum rules by 



for {jYf: 

for <k\mVf-- 



47^2 



(s ± 2ml) a/1 - 4m2/s 



27r2 



2mi^Jl — Aml/s 



:b.7) 



{ss) and to double the terms in ms{ss), m. 



In addition, one has to substitute 
and ms{sagGs) . 

We evaluate the sum rules using the input given in Table |Al The results are given in 
Sec. [23] and Tab. [H 



C Sum Rules for Twist-3 Matrix Elements 

The chiral-even twist-3 parameters C,\j^, , 1 ^Ik* ^^^l be determined from the correlation 
function 



afj, 



d^ye-^py{0\Tq{z)gGa,{vz)j,j,s{0)s{y)^,qiy)\0) = ipz)\2 - D)Ul^.{v,pz) ; (C.l) 



D is the number of dimensions. In terms of hadronic contributions, the correlation function 
is given by 



n 



3:K* 



{v,pz) 



m 



K* 



22 r 



(C.2) 
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the dots denote contributions from higher-mass states. The parameters Kg^*, (jjI^^* and Ag^^* 
can be obtained from an analogous correlation functions n|.^*(t>,p2;) with 

For the chiral-odd operator, one has 

i J d''ye-'Py{0\Tq{z)a,^gG,^ivz)s{0)siy)apMym = i(pzfui.j,4v,pz) . (C.3) 
All three correlation functions 11 can be written as 
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+ ^asTi5{ai)5{a2) {qq){ss). (C.6) 

Here we have dropped all terms that vanish upon Borelisation. The above expressions in- 
clude all terms to second order in the quark masses. One comment is in order concerning the 
contribution of the gluon condensate. Upon integration over a^, and subsequent expansion 
in powers of the quark masses, this contribution contains terms in ml ^ ln(m^ ^/ (— p^)), which 
are long-distance effects and must not appear in the short-distance operator product expan- 
sion of the correlation functions ( IC.ll) and ( 1C.3I) . As discussed in Ref. |34j, the appearance 
of these logarithmic terms is due to the fact that the above expressions are obtained using 
Wick's theorem to calculate the condensate contributions, which implies that the conden- 
sates are normal-ordered: (O) = (0| : O : |0). Recasting the operator product expansion in 
terms of non-normal-ordered operators, all infrared sensitive terms can be absorbed into the 
corresponding condensates. Indeed, using [M] 

2 

77? 77? ry 
{Q\-sgGs\Q) = (0| : -sgGs : |0) + In ^ (0| : ^ : |0), 

and the corresponding formula for q quarks, all terms in In tti^ ^ can be absorbed into the 
mixed quark-quark-gluon condensate and the resulting short-distance coefficients can be 
expanded in powers of m^^. In calculating the sum rules, we hence will use 

2 9 

In — > In — — . 

The QCD sum rules for the three hadronic parameters k^j^*, oj^^, and A^;^* describing 
the DA ^^-K*-, Eq. (13.111) . to NLO in conformal spin then read: 

{fk^y<^e-<*'^' I Xi^, = e-A^^ j Va {a, - ^) i I^s-^kM) ■ (C.7) 

The formulas for the other parameters are analogous. 

We evaluate the above sum rules in the Borel window = 1 GeV^ to 2.5 GeV^ and 
using the following values of the continuum threshold Sq: 

4(p) = (1.3±0.3)GeV^ = (l•3±0.3)GeV^ 4(</') = (1-4 ± 0.3) GeV^ 

4(p) = (1.5 ± 0.3) GeV^ , s^{K*) = (1.6 ± 0.3) GeV^ , s^{(f)) = (1.7 ± 0.3) GeV^ .(C.8) 

Numerical results, including the uncertainties from the variation of M^, sq and the input 
parameters of Tab. \^ are given in Tab. [TJ 
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D Loop Integrals 



For the interested reader, we collect the loop integrals needed for calculating the correlation 
functions in App. O 

At one loop, one has (with = 0) [15] : 



/ 



where the integration measure is defined as d^k = i/^AnY \d^k\ and fk is an arbitrary 
numerical factor. 

One also needs the following integral: 



T{c}T{d) Jo 
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2T{c)T{d) 



X / dMe*(i-")^'^'^M^/2-'^(l - (j + tfiil - u){kz)) . (D.2) 







Two-loop integrals are obtained by combining the above one-loop integrals. To obtain the 
correlation functions as integrals over T>a, one has to perform a variable transformation from 
{u,w) to (ai,a2), where the precise transformation relations are fixed by the "canonical" 
form of the exponential: exp{—ipz{a2 + v{l — ai — 0:2) })• Note that flD.2p contains an extra 



factor kz in the last line which comes from the Taylor expansion of the exponential. This 
factor can be made disappear by partial integration of the final result and hence does not 
appear in the explicit formulas for the correlation functions given in App. [Cl 
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